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ABSTRACT 

We present an observationally motivated model to connect the AGN and galaxy populations at 0.2 < 
z < 1.0 and predict the AGN X-ray luminosity function (XLF). We start with measurements of the 
stellar mass function of galaxies (from the Prism Multi-object Survey) and populate galaxies with 
AGNs using models for the probability of a galaxy hosting an AGN as a function of specific accretion 
rate (the rate of supermassive black hole growth scaled relative to the host stellar mass). Our model is 
based on measurements indicating that the specific accretion rate distribution is a universal function 
across a wide range of host stellar mass with slope 71 w —0.65 and an overall normalization that 
evolves strongly with redshift. We test several simple assumptions to extend this model to high 
specific accretion rates (beyond the measurements) and compare the predictions for the XLF with 
the observed data. We find good agreement with a model that allows for a break in the specific 
accretion rate distribution at a point corresponding to the Eddington limit, with a steep power- law 
tail to super-Eddington ratios with slope 72 = — 2.lj^'g. We convert between specific accretion rate 
and Eddington ratio by assuming a scaling between black hole mass and host stellar mass with an 
intrinsic scatter of ±0.38 dex. Our results show that samples of low luminosity AGNs are dominated 
by moderately massive galaxies {M.* ~ lO 10 " 11 ^©) growing with a wide range of accretion rates 
- a consequence of the shape of the galaxy stellar mass function rather than a preference for AGN 
activity at a particular stellar mass. The observed population of the most luminous AGN may be 
severely skewed to the most extreme sources with elevated black hole masses relative to their host 
galaxies and in rare phases of very rapid accretion. 
Keywords: galaxies: active - galaxies: evolution - X-rays: galaxies 



1. INTRODUCTION 

The universe has evolved rapidly over the last ~ 8 bil- 
lion years, since a redshift z ~ 1. The total star forma- 
tion rate density, the rate at which new s tars are being 
forme d, has dropped by a factor > 10 fe.g.|Wilkins et al.1 
120081 iZhu et all [211091 iRujopakarn et allboiOft . A sim- 
ilar decline is seen in the total rate of accretion onto 
supermassive black holes (SMBHs), which is tracked by 
the l uminosity density of Active Galactic Nuclei (AGNs) 
(e.g. Barger et al. 2005; Silv erman et al. 2008; Aird et all 
12010 ). This correlation suggests that the processes reg- 
ulating the formation of stars throughout galaxies are 
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somehow related to the processes that drive gas into their 
very central regions, feeding the growth of SMBHs and 
prompting periods of AGN activity. It is thus vital to un- 
derstand how the rapid decline of SMBH accretion since 
z ~ 1 is related to the co-evolving properties of their host 
galaxies. 

The AGN luminosity function provides the principal 
tracer of the distribution of SMBH accretion over the 
history of the universe. A variety of wavebands and iden- 
tification techniques have been used to measure the lu- 
mino sity function of AGNs out to high redshifts (z ~ 6) 
(e.g. IBovle et all 119871: IPage et al.l 119971: IRichards et al l 
2001 JAssef et al 120111: lUeda et all 120031: lEbrero et al.1 
2009; Air d et al.l I2010D . These studies have revealed a 



"downsizing" beha vior in the evolution of AGNs (e.g. 
iBarger etT al. 2005). This downsizing is characterized by 
a rapid decline in the number density of the most lumi- 
nous AGNs since a peak at z ~ 2, while the number den- 
sity of lower luminosity AGNs evolves much more weakly 
and peaks at lower z. 

To understand the underlying physical processes 
that drive the observed evolution of the AGN lu- 
minosity function it is vital to connect the AGN 
population to the galaxies that host them. A 
number of studies have indicated that AGNs are 
preferenti ally found in t he most massive galax- 
ies (e.g. iKauffmann et all 120031: iDunlop et alJ l2003t 



Schawinski et al.l 120071: iNandra et all 120071: iCoil et all 



20091: iHickox et al.1 120091 iBongiorno et all 1 2 01 2ft . These 
galaxies tend to have red optical colors and very low lev- 
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els of current star formation. Such trends suggest an 
association between the presence of an AGN and the 
quenching of star formation throughout a galaxy, pos- 
sibly due to feedback from the AGN. A number of more 
recent studies, however, have shown that AGN hosts have 
a similar distributio n of colors to galaxies of equivalent 
stella r masses (e.g. Kue et all l2010t iCardamone et al.l 
l201Cft. 

lAird et al.l (|2012f ) recently showed that the predom- 
inance of AGNs in more massive galaxies is in fact a 
selection effect; the probability of a galaxy hosting an 
AGN is determined by a power-law distribution of spe- 
cific accretion rates — the rate of accretion scaled by the 
stellar mass of the host — where the distribution itself 
does not depend on stellar mass. AGNs are present 
in galaxies with a wide range of stellar masses, but 
those in more massive gal axies are simply more lumi- 
nous. Ha inline et al.l (|2012h and Bongior no et al.l ()2012f ) 
have presented similar findings, extending the investiga- 
tions to higher redshifts. Nevertheless, these results indi- 
cate that the overall distribution of the stellar masses of 
galaxies, traced by their stellar mass function, is a vital 
component in determining the distribution of observed 
AGN luminosities, traced by their luminosity function. 

In this paper we develop a simple, observationally mo- 
tivated model that connects the evolution of AGNs and 
their host galaxies at 0.2 < z < 1.0. We start with the 
observed stellar mass function of galaxies and populate 
the galaxies with AGNs using a model for the probabil- 
ity of hosting an AGN as a function of specific accre- 
tion rate. Our m odel is based on measurements from 
lAird et all (|2012f) at low-to-moderate accretion rates. 
We extend the model to higher accretion rates using a 
number of well-motivated assumptions. We thus pre- 
dict the overall AGN luminosity function, which we then 
compare with the observed X-ray luminosity function of 
AGNs. This comparison allows us to test and refine our 
model, providing important insights into the nature of 
the AGN population and the physical underpinnings of 
the AGN luminosity function. We describe our observa- 
tional datasets in Section [2] and use these to construct 
our model in Section [3] Section |4] discusses our results 
and we state our overall conclusions in Section [5] 

2. OBSERVATIONAL DATASETS 

2.1. The galaxy stellar mass function 

To track the overall galaxy population we adopt recent 
meas urements of the g a laxy s tellar mass function (SMF) 
from iMoustakas et al.l (|2013t ). The SMF was measured 
using PRIMUS, a 9.1 deg 2 low-resoluti on spectroscopic 
redshift survey of galaxies to i ~ 23 (jCoil et al.l 120111 : 
iCool et al.i r2013D From the parent samp le of ~ 120, 000 
robust redshifts, IMoustaka s et al.1 ([2013) constructed a 
flux-limited sample of ~ 40, 000 galaxies at 0.2 < z < 1.0 
over an area coverage of 5.5 deg 2 . Stellar masses were de- 
termined using fits to the observed spectral energy dis- 
tributions (SEDs) based on ultra-violet, optical, near- 
infrared and mid-infrared photometry from GALEX, 
Spitzer /IR AC and a range of groun d-based imaging cam- 
paigns (see IMoustakas et al.ll2013l for full details). The 
SMF was constructed using the V max method, carefully 
accounting for the PRIMUS targeting weights and red- 
shift success rates, in six redshift bins between z = 0.2 



and z = 1.0 (corresponding to roughly equal intervals 
of cosmic time). For this study we adopt Schechter fits 
to the total SMF in each individual redshift bin, fully 
propagating the covariant uncertainties in the model fit 
parameters. A single Schechter form may not provide 
a full description of the shape of the SMF, but it does 
provide a good fit for the range of redshifts and stel- 
lar masses considered here. The Schechter fits allow us 
to smooth between the V m&K measurements in different 
stellar mass bins and perform moderate extrapolation to 
lower and higher stellar masse s. 

The Moustakas et al. (2013) results indicate remark- 
ably little evolution in the global SMF of galaxies in this 
redshift range: the cumulative number density of > 
1O 1O A10 galaxies has increased by ~ 30% since z ~ 0.6 
and the cumulative number density of Al* > 1O 11 A^0 
galaxies has changed by < ±10% since z ~ 1. Much 
stronger, differing evolution is found when the sample is 
divided into quiescent and star-forming galaxies, placing 
strong constraints on the rates of star formation quench- 
ing that builds up the quiescent population. In this pa- 
per, however, we consider only the global SMF and at- 
tempt to reconcile the relative lack of evolution in the 
SMF of galaxies with the strong evolution observed in 
the AGN population. 



2.2. The distribution of specific accretion rates 

To populate galaxies with AGNs we adopt measure- 
ments of the distribution of SMBH a ccretion rates for a 
given stellar mass and redshift from lAird et al.l (120121) . 
In lAird et al.l (|2012D we identified AGNs within the 
PRIMUS galaxy sample using hard (2-10 keV) X-ray 
observations of three of the PRIMUS fields. We then 
studied the probability of a galaxy hosting an AGN as a 
function of the stellar mass. 

We used hard X-ray selection to identify AGNs be- 
cause X-ray emission from a moderately luminous AGN 
(L2-iokcV 10 42 erg s _1 ) will dominate over other pro- 
cesses within the host galaxy (such as the combined 
emission from stellar X-ray binaries). In addition, hard 
X-rays (> 2 keV) are able to penetrate moderate col- 
umn densities of obscuring material at these redshifts 
(equivalent hydrogen column densities Njj < 10 23 cm -2 ). 
Thus, hard X-ray selection provides a reasonably accu- 
rate probe of the AGN luminosity and should track the 
bulk of the AGN population. The effect of the (highly 
variable) flux limits of the X-ray observations over the 
three PRIMUS fi elds was fully accounted for by the 
lAird et al.l (|2012ft work. We note that hard X -ray se- 
lection will fail to identify the most heavily obscured, 
Compton-thick sources. It remains unclear how large a 
fraction of the population is hidden by such high levels of 
obscuration. However, our measurements of the overall 
X-ray luminosity function (which we compare with our 
model predictions, see Section 12.31 below) also use hard 
X-ray selection to identify AGN and thus are subject to 
t he same inc o mplet eness . 

lAird et ahl (J2012) showed that the probability of a 
galaxy with a given stellar mass and redshift hosting an 
AGN could be described by a single power law distri- 
bution of specific accretion rates. This distribution is 
independent of stellar mass and evolves strongly with 
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Figure 1. Comparison of the binned measurements of the hard (2—10 keV) X-ray luminosity function (XLF) from Aird et al. ( 2010) (black 
circles, recalculat ed for the redshift bins shown here) and parametric model fits from various studies. The grey dashed line indicates the 
Aird ct al. (2010) LADE model evaluated at z = 0.25 (the center of the first redshift bin) and i s replicated in eve ry pan el. We also show the 
best fit Lu minosity-Dependent D ensity Evoluti on (LDDE) param eterisations of the XLF from Aird ct al. (2010, cyan) , [Ebrcro ct al, ( 120091 . 
magenta), Silverman ct al. (2008, orange) and Ueda et al. (2003, dark blue). While the details of the evolutionary behavior differ between 
the various parameterisations, all are approximately consistent with the observed data (black points). The evolution of the XLF from z 1 
to z ss 0.2 is predominantly driven by a shift of the XLF to lower luminosities with decreasing redshift (i.e. luminosity evolution). We 
adopt either the observed data points and their uncertainties (where available) or the range of the 5 different parameterisations to illustrate 
current measurements of the XLF and the uncertainties in subsequent plots. 



redshift, 

p(X I A4„, z) dlogA oc X 1 (1 + zf dlogA (1) 

where A oc Lboi/A'l* denotes the specific accretion 
rate, the rate at which mass is accreted by the SMBH 
(traced by the AGN bolometric luminosity, £ bol ) scaled 
by the total stellar mass of the host galaxy, and 
p(X | M.*, z) is the probability density per logarithmic in- 
terval in A of a galaxy at a given redshift hosting an AGN 
with specific accretion rate, A. We measured the slope of 
the probability distribution function as 7 = —0.65 ± 0.04 
and a strong redshift dependence, f3 — 3.5 ± 0.5. Our 
work spanned a wide range of stellar masses (9.5 < 
\ogM*/ Mq < 12) and focussed on low-to-moderate lu- 
minosity AGNs (42 < logix / erg s _1 < 44). Thus our 
results spanned a limited range in specific accretion rate, 
— 3 < log A < — 1, where A is defined relative to the Ed- 
dington limit assuming a single scaling between SMBH 
mass and the total stellar mass of the host galaxy, 
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These findings form the basis of our study in this paper. 
We investigate different models (see Section [3] below) to 
extend the A distribution to higher accretion rates and 



compare with the AGN luminosity function over a wide 
luminosity range. 

Broad-line QSOs were identi fied using the PR IMUS 
spectra and excluded from the lAird et all ((2012) sam- 
ples. In these sources the AGN emission dominates 
over the host galaxy at optical wavelengths, preclud- 
ing measurements of the stellar mass of the host galaxy . 
Over the X-ray luminosities probed bv lAird et al. (2012), 
Lx = 10 42 ~ 44 erg s _1 , the broad-line QSOs correspond 
to < 20% of the AGN population. Howev er, the broad- 
line fraction increases wit h luminosity (e.g. iSteffen et al.l 
120031 : iTreister et al.ll200"9T ) and thus excluding broad-line 
Q SOs may bias the slope of the A distribution measured 
by I Aird et al.l (|2012D . On the other hand, previous stud- 
ies have found that QSOs may be predominantly high 
accretion rate sources, accreting close to their Edding- 
ton limit (e. g. iKollmeier et al.N2006HTrump et al.ll201lL 
but see also Kelly fc Shenll2012| ) and dominate the pop- 
ulation at Lx > 10 44 ergs _1 . Our models (described 
below) that extrapolate the A distribution to higher ac- 
cretion rates may therefore account for the bulk of the 
broad-line QSO population. We discuss remaining issues 
regarding the exclusion of broad-line QSOs in Section HJ 

2.3. The X-ray luminosity function of AGNs 

The AGN luminosity function traces the overall distri- 
bution of accretion activity in AGNs in terms of the ob- 
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served luminosities at a particular waveband. Any model 
that populates galaxies with accreting SMBHs must be 
consistent with this observational constraint on the to- 
tal population. In this paper we compare our model to 
measurements of the 2-10 keV X-ray luminosity function 
(XLF) of AGNs at 0.2 < z < 1.0. 

In Figure H] we show binned estimates of the XLF from 
lAird et all ([2010D . recalculated in the redshift bins used 
for this study (black points with error bars). This work 
applied a sophisticated Bayesian methodology to account 
for uncertainties in photometric redshifts and X-ray flux 
estimates, track the X-ray sensitivity, correct for Edding- 
ton bias, and correct for incomplete optical identifica- 
tions. The binned estimates shown i n Figure [1] are calcu - 
lated using the N ^ s /N m( ii method (Mi vaii et al.ll200il ). 
which corrects for the effects of the X-ray sensitivity over 
the width of the bin. Error bars represent ler equiv- 
alent Poissonian errors for the number of objects in a 
bin. We also show a number of best-fit parametric mod- 
els, including both the luminosity-depe ndent densi t y evo - 
lution (LDDE) model determined by lAird et al.l (|2010f ) 
(cyan dashed line) and the conceptually simpler Lumi- 
nosity And Density Evolution (LADE) model where the 
XLF evolves in both luminosity and density but retains 
the sa me shape at all redshifts (green line). lAird et all 
(2010) concluded that it was not possible to statisti- 
cally distinguish between these models and thus advo- 
cated the simpler LADE model. In additio n, we show 
the b e st fit LDDE param eteriz ations from lUeda et al.1 
(2003). lEbrero et al.l (|2009h . and lSilverman et al.l (|2008). 

All of the model parameterizations demonstrate the 
same basic evolutionary trend, namely a reduction in the 
number density of more luminous AGNs, which can be 
characterized by an overall shift of the AGN population 
to lower luminosities as cosmic time progresses. This 
luminosity evolution dominates th e evolution descri bed 
by the LADE parameterization of lAird et al.l (|2010f) at 
these redshifts (z < 1) and has also been quantified by 
pure luminosity ev olution models in earlier works (e.g. 
iBarger et al.1 12005). However, the LDDE parameteriza- 
tions indicate that the shape of the XLF may be changing 
with redshift. Given these differences in the parameter- 
izations of the XLF, we choose to adopt the observed, 
binned data points and their uncertainties when com- 
parin g with our model predictions. Outside the range 
of the lAird et all (poToT ) data (L x > 10 45 erg s~ 4 at these 
redshifts) we conservatively take the minimum and maxi- 
mum of the 5 model parameterizations shown in Figured] 
to represent current uncertainties in the XLF. This XLF 
estimate is shown by the light-green hatched region in 
subsequent plots. 



3. AN OBSERVATIONALLY MOTIVATED MODEL TO 
CONNECT THE AGN AND GALAXY POPULATIONS 

In this section we present a simple, observationally mo- 
tivated model to populate galaxies with AGNs. We start 
with our measurements of the galaxy SMF and convolve 
with a model for the probability of hosting an AGN as 
a function of specific accretion rate. This convolution 
provides us with a prediction for the XLF of the AGN 
population, which we can compare with the observations. 



Our model for the XLF can be expressed as 

ip(L x ,z) = <f>(M*,z) *p(X | M*,z) (3) 

4>(M*,z) p(A(L x , M*) | M*,z) d\ogM* 



where 4»(M*,z) is the galaxy SMF and p(X \ M*,z) 
describes the probability density per logarithmic A in- 
terval for a galaxy of given AA* and z hosting an 
AGN with a specific accretion rate, A. This equa- 
tion allows us to predict the XLF, ip(Lx,z), assum- 
ing the luminosity-d ependent bolometric corrections of 
Hopk ins" et al.l (|2007t ) to convert from A to an X-ray lu- 
minosity, Lx, for a given M.*. In Figure [2] we show our 
three models for p(X | A4*,z), all of which are based 
on th e power-law distribution measured by lAird et aLl 
<|2012f ) that is independent of stellar mass but evolves 
strongly with redshift (see Equation Q] and Section 12.21 
above) . 

3.1. Model A 

We begin with a simple model (A), where we assume 
a single scaling applies between SMBH mass and the to- 
tal stellar mass of the host galaxy at all redshifts. We 
take the scaling between SMBH mass and galax y mass 
as M hh » 0.002A1* from lMarconi fc Hunt! (|2003D . where 
we have further assumed that 7W* ~ A4 bulge- Thus, we 
can directly translate the specific accretion rate into an 
Eddington ratio for the SMBH (see Equation [2]) . With 
this scaling assumed, we apply a strict cut to our spe- 
cific accretion rate distribution at the Eddington limit, 
log A = 0, as shown in Figure [5] (left panel). The grey re- 
gion tracks the uncertainty in our specific accretion rate 
distribution. We estimate the uncertainties using the 
errors in the best-fit model parameters from lAird et al.l 
(|2012| ). However, these errors assume a fixed model form 
and do not fully account for the observational uncertain- 
ties in the measurements of the A distribution at a given 
stellar mass or redshift. We therefore increase the un- 
certainty in the overall normalization by ~ ±0.2 dex, 
which provides better agreement between our model un- 
certainty and the errors in the indiv idual binned mea - 
surements of p(X \ M*,z) in figure 7 of lAird et al.l (2012). 
This increase also allows for additional systematic uncer- 
tainties in our model, such as the possible contribution of 
unobscured sources. In Figure [3] we compare the predic- 
tion assuming our model A for the specific accretion rate 
distribution (thick black line with the solid grey region 
indicating the uncertainty) with the observational con- 
straints on the XLF (green hatched region). Our simple, 
observationally motivated model is able to successfully 
reproduce the faint-end of the XLF at all redshifts to 
z ~ 0.7. The predominant evolution with redshift is 
accounted for by the strong evolution in the normaliza- 
tion of p(X | M*,z), which produces a change in the 
normalization of our predicted XLF. However, we fail 
to reproduce the correct behavior at the bright end of 
the XLF (Lx > 10 44 5 erg s" 1 ). Our model predicts an 
exponential decline in the XLF at high luminosities, a 
consequence of the Schechter form of the galaxy SMF 
and the sharp cutoff we apply in the specific accretion 
rate distribution. We note that our model does not re- 
produce any of the (relatively minor) changes that may 
be observed in the evolution of the faint-end slope of the 
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Figure 2. Assumed distributions of specific accretion rates (for a sample of galaxies at a given stellar mass and redshift), based on 
measurements (at -3 < log A < -1) bv lAird et"aT] H2012I1 . The units assume A = L bol (l.3 X 10 38 erg s _1 X 0.002-£i^ . In model A we 

assume a single scaling between the SMBH mass and stellar mass of the host galaxy; thus, A directly translates into the Eddington ratio 
in the assumed units and we apply a strict cutoff at the Eddington limit, log A = 0. The lines indicate the distribution at z = 1.0 (dotted), 
z = 0.6 (solid) a nd z = 0.2 (dashed) showing the evolution in the probability of a gala xy hosting an AGN with redshift as measured by 
Aird ct al. (2012). The grey region indicates the uncertainty in the model fit to the Aird ct al. (2012) measurements, including an additional 
±50% uncertainty in the normalization. In model B we allow for an intrinsic lognormal scatter of 0.38 dex in the scaling between SMBH 
mass and total host stellar mass. This scatter softens the cutoff in terms of specific accretion rate and predicts a population of SMBHs 
with high masses relative to their host accreting close to their Eddington limit, which thus have high specific accretion rates. In model C 
we also allow for a steep power-law distribution to higher specific accretion rates, corresponding to a rare population of SMBHs accreting 
at rates above their Eddington limit. 

XLF with redshift nor any change in the break of the 
XLF. At z > 0.7 we underpredict the number density of 
L x « io 4 3-5-44.5 erg s -i AGNs, and the observed faint- 
end slope is somewhat flatter than our prediction. 

3.2. Model B 

In our second model (B), we adopt a more realistic 
scaling relation between the SMBH mass, TWbh, and 
host stellar mass, that allows for intrinsic scat- 

ter in the relationship. We adopt an intrinsic scatter 
of 0.38 dex, consistent with measurements of the local 
relationship b etween sphero i d stella r mass, Al spheroid; 
and A^bh by iBennert et al.1 (|2011a[ ). This work used 



updat ed dynamical A'lbh estimates from iGultekin et al.1 
(2009) for a sample of 18 local elliptical and SO galaxies; 
spheroid stellar masses were es timated using the J, H 
and K magnitudes measur ed by iMarconi" fc Hunti |2003) 
and the lAuger et al.l (|2009t ) Bayesian stellar mass es tima- 
tion code. We note that the IBennert et all (|2011aj ) rela- 
tionship was measured using the spheroid (bulge) stellar 
masses, but we assume that the same intrinsic scatter 
can be applied to the total stellar mass - SMBH mass 
relationship. 

In model B we assume this 0.38 dex scatter in the AJbh- 
M* relationship to convert between A and Eddington ra- 
tio, applying a sharp cutoff in the probability of a galaxy 
hosting an AGN above the Eddington limit. The scatter 
in A'lbh-A'f* results in a much softer cutoff in terms of 
specific accretion rate, A (see Figure center). The tail 
in p(X | A4», z) at high specific accretion rates (log A > 0) 
can be attributed to the probability of finding galaxies 
with elevated SMBH masses relative to the host (due to 
the intrinsic scatter) accreting close to their Eddington 
limits. 

Figure [4] compares the observed XLF and our predic- 
tion using model B. The softer cutoff in the A distri- 



bution due to the scatter of the A^bh-AI* relationship 
results in a smoother break in our predicted XLF com- 
pared to model A (Figure [3]). Furthermore, at the bright 
end of the XLF (L x > 10 44 5 erg s" 1 ), AGNs with typ- 
ical mass SMBHs are outnumbered by the small frac- 
tion of lower stellar mass galaxies with elevated SMBH 
masses, due to the intrinsic scatter in the A^bh-Al* re- 
lation and the steep slope of t he SMF at high stellar 
masses (M* > 10 n A4 o , see also lSomerville] |2009'). This 
effect leads to a somewhat flatter predicted bright-end 
slope than found with Model A. Nonetheless, our predic- 
tion still declines more rapidly than the observed XLF at 
Lx > 10 45 erg s _1 , especially at high redshifts (z > 0.5). 

3.3. Model C 

Prompted by the lack of agreement at the bright end 
of the XLF using model B, our third model (C) re- 
moves the requirement of a sharp cutoff at the Edding- 
ton limit and instead adopts a steep power-law distri- 
bution with a tail to very high Eddington ratios, while 
retaining the 0.38 dex scatter in the Mbh~-M* scaling. 
The Eddington limit only strictly applies to a spheri- 
cal geometry and studies of accretion disk physics do al- 
low f or scenarios when the Eddington limit is violated 
(e.g. Uaroszvnski et~aLl 119801: lAbramowicz et al.l 119881 : 
iBegelman et al.l 120061 ). Indeed, SMBHs accreting above 
their Eddington limits are rare but are found in suffi- 
cient l y large samples of A GNs (e.g. iMcLure fc Dunlopl 
l20f)l IKellv fc Shed 12011 . The slope of this high ac- 
cretion rate tail has not been directly measured in our 
observed specific accretion rate distri bution (although 
hints of such a turnover may be seen by iBongiorno et all 
2012). We therefore introduce this as a single free pa- 
rameter, 72, in our model that we fit using the original 
binned data points in our observed XLF (shown in Figure 
[T]). We constrain the slope to be 72 = — 2.lto's- 
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log L x / erg s" 1 

Figure 3. Predictions for the XLF (black solid line, grey region shows propagated lcr uncertainty in the model) at a range of redshifts 
based on the galaxy SMF and assuming model A for the distribution of accretion rates. Model A assumes the mass of the SMBH has a 
single scaling with the galaxy stellar mass and thus the distribution of specific accretion rates has a power-law form with a sharp cutoff 
corresponding to the Eddington limit. The colored lines show the contributions to the XLF from galaxies over limited ranges of stellar 
mass. We compare to direct observations of the XLF (green hatched regions, which track the uncertainty in the observations — see Figure 
[T] and Section 12.31 for details). Our model successfully predicts the faint end of the XLF (Lx $ 10 44 erg s — 1 ) but falls over much more 
rapidly than the observed slope for more luminous AGNs and fails to reproduce the evolution in the break of the XLF. 



In Figure [5] we compare our prediction from model C 
with the observed XLF. With model C we are able to 
produce an XLF that is in good agreement with the ob- 
served data (considering the uncertainties in both our 
prediction and the observations). However, there are 
minor discrepancies between our observations and our 
model for the XLF. Our model underpredicts the num- 
ber density of L x w io 43 - 5 - 44 - 5 er g s" 1 AGNs at z > 0.7, 
although the discrepancy is only at the ~ 2er level when 
considering the uncertainties in both the observed XLF 
and our model prediction (propagated from the SMF 
and measurements of the specific accretion rate distri- 
bution). In addition, the evolutionary behavior of our 
model is predominantly a density evolution. We do not 
reproduce the luminosity-dependent evolution found in 
a number of studies of the XLF, although our model 
prediction does provide reasonable agreement with the 
observed XLF considering the current uncertainties. We 
discuss these issues in Section [4] 

Our observationally motivated model allows us to in- 
terpret the observed AGN population at 0.2 < z < 1.0, 
as traced by the XLF. As shown by the colored lines in 
Figure El our model predicts that the low Lx popula- 
tion is dominated by AGNs in hosts with stellar masses 
~ 10 10 - 5-11 A^o accreting over a wide range of specific 
accretion rates (consistent with the predomina nce of such 
host galaxies in X-ray selec ted samples, e.g. Kue et afl 
l2010t [Bongiorno et al J 120121 ). Around the break in the 
XLF the population is dominated by moderately massive 



galaxies ~ 10 11 A / J©) accreting close to their Ed- 

dington limit. The bright end of the XLF is dominated 
by massive galaxies (Ai* > lO 11 A / (0) accreting close to 
or above the Eddington limit. Assuming a 0.38 dex scat- 
ter in the intrinsic AJbh — M* relation, we predict that 
80% of detected AGN with L x > 10 44 ' 5 erg s" 1 will have 
SMBH masses that are above the fiducial Mhh — M* re- 
lation; the predicted median A4bh/ 'M* for this luminous 
population is ~ 0.25 dex higher than the median for the 
overall galaxy population. This effect can be explained 
by the steep slope of the SMF at high masses and the 
intrinsic scatter in the M.bh-M-* relation. It does not 
indicate that SMBH masses are intrinsically higher in 
AGN host galaxies. Finally, the evolution in the normal- 
ization of the XLF is explained by a rapid drop in the 
probability of a galaxy hosting an AGN at later times for 
all specific accretion rates and across all stellar masses. 

4. DISCUSSION 

We have shown how with a simple, observationally mo- 
tivated model we are able to explain the predominant fea- 
tures of the XLF of AGNs and its evolution over 0.2 < 
z < 1.0. We start with the observed SM F of galaxies (and 
its ve ry mild evolution) from PRIMUS (Moustaka s et all 
2013) and populate these galaxies with AGNs based on 
a model for the probability of hosting an AGN as a func- 
tion of specific accretion rate. Our model is based on the 
observations at low-to-moderate accretion rates (equiva- 
lent to Eddington ratios ~ 10~ 3 — 0.1) and over a wide 
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41 42 43 44 45 46 42 43 44 45 46 42 43 44 45 46 

log L x / erg s" 1 

Figure 4. As Figure [3] but assuming model B for the distribution of specific accretion rates. We account for intrinsic scatter of 0.38 dex 
in the A4bh — relation to calculate the Eddington limit, where we truncate the probability of a galaxy hosting an AGN. This scatter 
softens the cut off in the distribution of specific accretion rates and flattens the predicted bright-end slope of the XLF. Our prediction still 
falls off more rapidly than the observed XLF at high Lx- 




Figure 5. As Figure[4] but assuming model C for the distribution of specific accretion rates. In addition to the scatter in the Mbh — A4* 
relation, in this model we also include a steep power law tail to high values of A, corresponding to super-Eddington accretion rates. Our 
model is able to reproduce the dominant features of the observed XLF from z ~ 1 to z ~ 0.2. 
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range of stellar m asses (9.5 < log M*/Mq < 12) from 
lAird et all (|2012l ). We extend our model to higher accre- 
tion rates (equivalent to Eddington ratios > 0.1) using a 
small number of straightforward assumptions. Our final 
model has one free parameter — the slope of the distri- 
bution at super-Eddington accretion rates — which we fit 
using the observed bright end of the XLF. These results 
provide new important insights into the nature of the 
AGN population and its evolution since z ~ 1. 

4.1. The nature of the AGN population 

At low X-ray luminosities (the faint end of the XLF), 
the number density of AGNs is dominated by moderately 
massive galaxies (A4» ~ lO 10 5-11 ^©) with a range 
of specific accretion rates given by the pow er-law dis- 
tribut ion with slope w —0.65 measured by lAird et al.l 
<|2012f ) . The predominant contribution by moderate mass 
galaxies is due t o the shape of the SMF. As shown in 
lAird et all (|2012f ). AGN activity does not appear to have 
a preference for a particular stellar mass range. In our 
model, galaxies with lower or higher stellar masses are 
just as likely to host an AGN with the same distribution 
of accretion rates, but these sources do not constitute a 
significant fraction of the overall AGN number density. 

The bright end of the XLF probes a different regime 
of the AGN population. In our model C, we propose a 
break in the accretion rate distribution at a point corre- 
sponding to the Eddington limit and a steep power-law 
tail that allows for a rare population of sources with very 
high (super-Eddington) accretion rates. We also allow 
for an intrinsic scatter in the scaling between the SMBH 
mass and the stellar mass of the host galaxy. We find 
that the bright end of the XLF consists mostly of mas- 
sive galaxies (M* > 10 11 A / fo) accreting close to or above 
their Eddington limits. This population is dominated by 
AGNs with SMBH masses that are elevated relative to 
the fiducial Mbh-M* relation, which is a result of scat- 
ter in the scaling rel ation and th e steep decline of the 
SMF at high masses. ILauer et al.l (|2007n discussed how 
such a bias could affect studies of SMBH mass scaling re- 
lations using broad line QSOs outside the local universe. 
In fact, recent studies have found tentative evidence for 
strong evolution in the scaling between .Mbh and host 
spheroid mass (e.g iTreu et al.N2007l: IJahnke et al.l l2"00l 
IMerloni et all 120101: IBennert et al.l 1201 lbf) . Our model, 
which assumes a constant scaling and scatter in the A4bh~ 
relation, may provide a simple explanation for such 
findings, although a full investigation is beyond the scope 
of this paper. 

We note that i n our study of sp ecific accretion rates 
(|Aird et al.ll20ll see Section 12.21) we did not measure 
p(X | J\A* : z) for high accretion rate sources (log A > 
— 1). In fact, our study excluded broad-line QSOs 
and was limited to moderate luminosity AGNs (Lx = 
10 42-44 erg s -i) We demonstrated in lAird et all (|2012h 
that the exclusion of QSOs did not significantly bias 
our measurement of the slope of the specific accretion 
rate distribution at —3 < log A < —1. In fact, QSOs 
have been shown to dominate the popu lation of lu- 
mmous AGNs (L x > 1 44 ergs" 1 , e.g. iBarger et all 
120051 : iTrump et al.l l2009h . Nevertheless, measurements 
of p(X | z) at log A > — 1, including QSOs and lumi- 
nous (Lx > 10 44 erg s _1 ) AGN, are vital to confirm our 



interpretation of the XLF and our model relating AGNs 
and their host galaxies. Such studies, however, present 
a number of difficulties due to the relative rarity of such 
sources, the strong biases that can affect the properties 
of the observed population, and the challenge of disen- 
tangling the properties of luminous AGNs and their host 
galaxies. 

4.2. The evolution of AGNs 

We attribute the evolution in the XLF since z ~ 1 to 
a reduction in the overall probability of a galaxy hosting 
an AGN as cosmic time progresses, while the shape of 
the distribution of accretion rates remains roughly con- 
stant. This behavior can be interpreted as a rapid drop 
in the rate of AGN triggering (the duty cycle) or a shift to 
lower characteristic accretion rates, either of which may 
be driven by a reduction in the availability of cold gas 
with cosmic time. Despite this rapid evolution, AGNs 
have a wide range of specific accretion rates (with a con- 
sistent power-law distribution) and are hosted by galax- 
ies across a wide range of stellar masses at any epoch. 

The details of the physical processes that trigger AGNs 
and give rise to this distribution of specific accretion rates 
are still unclear. The Eddington limit appears to set the 
maximum rate of SMBH growth at all times (although 
in rare or short-lived phases this limit may be violated). 
Furthermore, the power-law distribution of accretion 
rates may reflect variability in the level of AGN activity 
over relatively short timescales (compared to the lifetime 
of a galaxy), which could also be related to the physical 
processes in very central regions of t he galaxy that con - 
trol the rate of SMBH accretion fe.g. lNovak et al.N201lh . 
However, variations in larger-scale fueling mechanisms — 
such as mergers, disk instabilities, cosmological accre- 
tion of cold gas, or mas s loss from the stellar popu- 
lation of the galaxy (e.g. Kauff mann fc Haehneltl 120001 : 
IHopkins fc Hernquistl 120061; ICiotti fc Ostrikerl 12007ft — 
may also have an impact on the accretion rate distribu- 
tion. Indeed, changes in large-scale fueling mechanisms 
may drive the longer-term evolution that leads to lower 
characteristic accretion rates or reduced duty cycles at 
later cosmic times. The origin of this power-law distri- 
bution is is thus a key problem for theoretical studies of 
AGN physics. 

Our model is successful at predicting the predomi- 
nant evolution of the XLF. However, we underpredict 
the number density of moderate luminosity AGNs (Lx ~ 
10 43-44 erg g-i) at z > 7j although the discrepancy is 
only at the ~ 2a level. Furthermore, our model lacks 
the strong, luminosity-dependent evolution favored by 
most studies of the XLF over wid e r redshift ranges (out 
to z ~ 3, e.g. lUeda et all 120031: jHasinger et al.l 120051 : 
lEbrero et all 120091: lAird et al.ll2010D . These luminositv- 
dependent schemes are associated with "downsizing" of 
AGN activity. Observationally, downsizing is seen as a 
rapid evolution in the number density of luminous AGNs 
that peaks at z ~ 2, while the space density of lower lu- 
minosity AGNs evolves more weakly and peaks at lower 
z. Downsizing may also be reflected in a flattening of the 
faint-end slope of the XLF, possibly seen in t he observed 
XLF in our highest redshift bin (although lAird et all 
120101 concluded that a constant shape of the XLF was 
consistent with the observational data, carefully consid- 
ering the uncertainties). This downsizing behavior is 
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thought to reflect the build up of SMBH masses, whereby 
the highest mass SMBHs form and grow early in the his- 
tory of the universe, whereas lower mass SMBHs carry 
out most of their growth at later times. Evidence for this 
behavior may be seen in the growth rates of AGN s in the 
local universe (z < 0.3, e.g. lHeckman et aLl l2004'). 

Given the success of our model at reproducing the pre- 
dominant evolution of the XLF at 0.2 < z < 1.0, it ap- 
pears that any downsizing behavior may be a secondary 
effect at these redshifts. Our results indicate that the 
predominant evolution of the AGN population at z < 1 
can be associated with a rapid decline in the probability 
of hosting an AGN for galaxie s of any st e llar m ass as 
cosmic time progresses. In fact. lAird et al.l (|2012fi found 
that the distribution of specific accretion rates was con- 
sistent over a wide range of stellar masses — a key assump- 
tion of our model — and found no evidence of a downsiz- 
ing behavior in terms of host stellar mass. Nevertheless, 
the discrepancies between our model and the observed 
XLF suggest that the specific accretion rate distribution 
may evolve in a mass-dependent way. Further study is 
required to confirm such behavior and motivate any re- 
finements to our model. 

4.3. The relationship between AGN activity and star 
formation 

The evolution of the AGN population has strong par- 
allels with the star formation history of the universe, 
which may also undergo a downsizing behavior where 
the most massive galaxies a ppear to form earlier in the 
history of the universe (see iFontanot et "ail [2009. for a 
critical overview of the observational evidence). How- 
ever, it is now becoming clear that the rapid decline in 
the global star formation rate density since z ~ 1 — 2 is 
primarily due to a drop in specific star formation rates 
for all galaxies, with any stellar-ma ss-dependent down- 
sizing being a secondary effect (e.g. iNoeske et al.l 120071 : 
iZheng et a l. 2007; Karim et al.ll201lD . similar to our find- 
ings for AGNs. Nevertheless, how and why these pro- 
cesses are connected remains a key open question. 

Our simple observationally motivated model assumes 
that the probability of a galaxy hosting an AGN is deter- 
mined by a universal specific accretion rate distribution 
that is independent of host stellar mass or star formation 
properties. Thus, the observed AGN population (traced 
by the XLF) is determined by the galaxy SMF, the scal- 
ing between M.bh and Al*, and the redshift evolution of 
the AGN duty cycle. However, studies of the correlation 
between star formation and AGN accretion rates paint 
a more complex picture, indicating that star formation 
and AGN accretion rates ar e connected, at least for the 
most luminous sources (e.g. iNetzerl [20091 : iRosario et al.1 
120121) , when considering the nuclear regions of galaxies 
(e.g. iDiamond-Stanic fc Riekd I2012T ) or when averaged 
over a large populat ion of galaxies (e .g. iMullanev et al.l 
12012ft . Furthermore. lAird et al.1 (|2012[) did find evidence 
for a weak (factor ~ 2) enhancement in the probability 
of hosting an AGN for blue galaxies (indicating current 
star formation) compar ed to red (quiescent) galaxies, and 
iMoustakas et al.l (|2013[ ) found that the SMFs of the star- 
forming and quiescent galaxy populations follow substan- 
tially different evolutionary paths. Further study is thus 
required to fully unveil the interplay between AGN and 
star formation activity. 



4.4. Other models of AGN -galaxy co-evolution 

A number of prior studies have presented models to 
link the evolution of AGNs and galaxies, either from an 
observationally motivated standpo int (e.g. Shanka r et al.l 
120091 120131 : IConrov fc Whitdl2013l ) or incorporating pre- 
scriptions for SMBH growth into either l arge-scale hydro- 
dynamic cosmological sim ulations (e.g. iBooth fc Schayd 
120091 : iDeGraf et al.l 12012ft or semi-analytic galaxy for- 
mation models (e.g. iBower et al.l 120061 : iSomerville et al.1 
2008; Fanidakis_et^|2012|). Given the current uncer- 
tainties regarding SMBH fueling mechanisms and the 
physics of AGN accretion, the luminosity function is of- 
ten used to constrain the key parameters in the prescrip- 
tions for SMBH growth, such as the AGN duty cycle 
or the relati ve rates of different AGN triggering mech- 
amsms (e.g. iHopkins et all 120061: iShankar et all 120131 : 
IConrov fc White! 120131: IDraper fc Ballantvnd 120121 . 
However, up-to-date semi-analytic models, including so- 
phisticated prescriptions for SMBHs, are able to pre- 
dict t he AGN luminosity func tion (e.g. iFanidakis et al.1 
2012; Hirschm ann et al.1 12012T ). although it can be dif- 
ficult to assess the level of tuning required for a suc- 
cessful model or robustly test t he underlying physica l 
assu mptions. For example, t he IFanidakis et al.l (2012) 
and iHirschmann et all (|2012ft models both reproduce 
the observed AGN luminosity functions but infer very 
different levels of AGN triggering due to disk insta- 
bilities versus merger-dr iven inflows. Furthermore, the 
IHirschmann et all (|2012f) model reproduces the observed 
downsizing trends in the XLF by modifying their fiducial 
model to include a limit to accretion rates that depends 
on the cold gas fraction, as well as a prescription for AGN 
triggering due to disk instabilitie s and a "heavy seeding " 
prescription for the first SMBHs. Fanidakis et jil] (|2012ft . 
however, attribute the observed downsizing to obscura- 
tion effects, incorporated into the ir model using a pre- 
scription based on observations by iHasingeil (|2008|) . 

Our simple, observationally motivated model takes a 
very different approach. We adopt the observed SMF 
of galaxies and populate the galaxies with AGNs using 
models for p(X \ M*,z), based on actual measurements, 
requiring only a few straightforward assumptions to ex- 
tend to higher accretion rates. Our final model has one 
free parameter which is tuned to match the observed 
bright end of the XLF but does not alter our model 
predictions at lower luminosities. Unlike the models 
described above, our observationally motivated scheme 
does not provide a complete picture of the relationship 
between AGNs and their host galaxies or the underlying 
physical processes that control their co-evolution. How- 
ever, our results indicate that the universal distribution 
of specific accretion rates, independent of host stellar 
mass, is the key underlying property of the AGN popu- 
lation that is needed to reproduce the XLF. Thus, the 
distribution of specific accretion rates is a vital obser- 
vational constraint for theoretical models of AGNs and 
their evolution within a cosmological context. Our model 
also provides a framework to interpret different measure- 
ments, understand underlying observational biases that 
may skew the properties of samples of AGNs, and moti- 
vate future observations. In fact, the lack of downsizing 
in our model predictions and the discrepancy with the 
observed XLF for L x « 10 43 ~ 44 erg s" 1 at z ~ 0.7 - 1.0 



10 



AlRD ET AL. 



provides an important motivation for precise measure- 
ments of p(X | M.*, z) at these redshifts and higher. Ul- 
timately, conclusive tests of cosmological AGN-galaxy 
models may require comparisons of predictions in the 
multivariate space of galaxy stellar mass, star formation 
properties, and the distribution of AGN accretion rates. 

5. CONCLUSIONS 

We have developed a simple, observationally moti- 
vated model to link the evolution of AGNs and their 
host galaxies that predicts the X-ray luminosity func- 
tion of AGNs at 0.2 < z < 1.0. We use new measure- 
ments of the galaxy ste llar mass function from PRIMUS 
(jMoustakas et al] 12013) combined with a simple model 
to populate these galaxies with AGNs. Our model is 
based on our measurements showing that the probabil- 
ity of a galaxy hosting an AGN is defined by a power-law 
function of specific accretion rate (with slope 71 ~ 0.65) 
that depends only on redshift a nd is independent of host 
stellar mass (jAird et al.1 120121 ). We extend our model 
to higher X-ray luminosities and accretion rates by al- 
lowing for a break in the specific accretion rate function 
at a point corresponding to the Eddington limit. We 
also allow for scatter in the scaling between SMBH mass 
and the stellar mass of the host galaxy. To match the 
bright end of the XLF, we introduce one free parameter 
corresponding to the slope of the power-law tail of the 
distribution of accretion rates above the Eddington limit. 
We constrain this power-law slope to be 72 = — 2.1+q| 
by fitting to the observed XLF at z = 0.2 — 1.0. 

Our model provides a simple picture of the connection 
between AGNs and t heir h ost galaxies. Based on the 
results of Aird et al.l ()2012l ). AGN activity may be trig- 
gered in galaxies of any stellar mass, with the growth 
rate of the SMBH being defined by a single distribution 
of specific accretion rates. The observed rapid evolution 
of the AGN population since z w 1 is attributed to an 
overall reduction in the probability of any galaxy host- 
ing an AGN, indicating a severe reduction in the rate of 
AGN fueling as cosmic time progresses. While the phys- 
ical processes driving this behavior remain unclear, they 
appear to operate in a uniform manner over a wide range 
of stellar masses. Given our success at reproducing the 
XLF at z < 1 with this simple model, it appears that the 
universal distribution of specific accretion rates is the key 
observational constraint for theoretical models of AGN 
evolution. 

We find that the observed XLF is dominated in num- 
ber density by moderately massive galaxies « 
lO 11 ^©), which is due to the combination of the shape 
of the galaxy SMF and the shape of the specific accretion 
rate distribution. This predominance of moderately mas- 
sive host galaxies in X-ray selected AGN samples does 
not represent a preference for AGN activity at any par- 
ticular stellar mass range. The observed population of 
the most luminous AGN may be severely skewed to the 
most extreme sources with elevated SMBH masses rel- 
ative to their host galaxies and in short-lived phases of 
very rapid SMBH accretion. These sources are likely to 
provide a biased view of the connection between AGNs 
and the properties of their host galaxies. 

While our simple model reproduces the predominant 
evolution of AGN population since z ~ 1, we do not pre- 
dict any downsizing in SMBH growth with cosmic time. 



We do not fully reproduce the luminosity-dependent evo- 
lution in the space density of AGNs or possible changes 
in the faint-end slope of the XLF at z <~ 0.7 — 1.0. Thus, 
the specific accretion rate distribution may evolve in a 
mass-dependent manner, although any downsizing be- 
havior appears to be a secondary effect at z < 1. Our 
model also lacks known connections between the star 
formation in galaxies and the growth of their SMBHs. 
These discrepancies motivate not only the need for high 
precision measurements of the XLF over a wide range of 
redshifts but also further detailed observational studies 
of the distribution of AGN accretion rates within well- 
defined samples of galaxies divided by stellar mass, red- 
shift, and star formation properties. 
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